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Abstract 

In addition to well-motivated scenarios like supersymmetric particles, the so-called exotic matter 
(quirky matter, hidden valley models, etc.) can show up at the LHC and ILC, by exploring the 
spectroscopy of high mass levels and decay rates. In this paper we use QCD-inspired potential 
models, though without resorting to any particular one, to calculate level spacings of bound states 
and decay rates of the aforementioned exotic matter in order to design discovery strategies. We 
mainly focus on quirky matter, but our conclusions can be extended to other similar scenarios. 
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I. INTRODUCTION 



Since the beginning of accelerator physics, mass spectroscopy has been playing a leading 
role in the discovery of particle and resonance states, and understanding of the fundamental 
interactions in the Standard Model (SM). For example, the first signals of charm and bottom 
quarks were in fact detected through the formation of J/ip and T bound states. 

On the other hand, current colliders like the LHC, or the ILC in a farther future, will 
likely continue this discovery program beyond the SM. It is conceivable that new (super) 
heavy bound states can be formed and, contrary to e.g. the toponium system, their basic 
constituents are prevented from decaying before the binding is effective. The goal of this 
paper is to perfom a prospective study of the spectroscopy of such exotic massive states, 
by making several reasonable assumptions about the interacting potential among the new- 
physics constituents which may differ from standard QCD. Furthermore, we will estimate 
leptonic decay widths of very heavy bound states by making specific assumptions on the 
quantum numbers of constituents, although not in a comprehensive way. 

A. Exotic scenarios 



During the last years, minimal extensions of the SM containing additional heavy particles 
charged under a new unbroken non-Abelian gauge group G v with fermions Q,Q have been 
proposed under the general name of "hidden valley" models which is a very general 
scenario containing such heavy particles but as well new sectors of lightweight particles to 
be observed. In these models all SM particles would be neutral under such the new G v 
group, while new particles charged under G v but neutral under the SM group would show 
up if large energy scales are probed. Higher dimension operators, induced e.g. by a Z' or a 
loop involving heavy particles carrying both Gsm and G v charges, should connect both SM 
and new physics sectors through rather weak interactions. 

In particular, if the G v group corresponds to SU(3), the fermions in the fundamental 
representation have been recently named as "quirks" [2j or iquarks j^]. 4] Actually this 
theory can be viewed as a certain limit of QCD where light quarks are removed and the 
typical scale A where the new interaction becomes strong is much smaller than the heavy 
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flavor masses. More generally, such kind of scenario can be put in correspondence with a 
class of "hidden valley" models as pointed out in [1]. In these models, quirks are defined to 
be new massive fermions transforming under both G$m and a general ( not only SU(3) ) 
non-abelian gauge group G v . 

It has been considered so far in the literature that A is smaller than Aqcd, but as well 
the particular case of hidden valleys with quirks in which A is greater than Aqcd has been 
studied in the literature The name of infracolour is used to design of the new 

gluonic degrees of freedom when A is much smaller than the quirk mass; in this work we 
will refer to it as i-colour (i-QCD) in correspondence with the name quirk. 

B. Exotic long-lived bound states 



In this section we will focus on the bound states of quirks. The phenomenology of such 



bound states has first studied in |2 



, and later an analysis on the spectroscopy of these 



systems was done by the authors of 8|. 



It is well-known that the large mass of the top quark in the SM prevents toponium to be 
formed since the constituents quarks would decay away too fast. A criterion for existence 
of such bound states is that the binding energy should be larger than the total decay 
width. For heavy onia states beyond the SM, however, the situation could be different. 
In particular, in the case of quirks Q , |l3 1 , js)] , its decay is prevented from the conservation 
of a quantum number. Regarding the dynamics of these quirks, according to (2|, one can 
distinguish among three possible energy scales: the first one is the 100 eV < A < keV 
range, where the quirk strings are macroscopic; the mesoscopic strings can be find at the 
range keV < A < MeV (which is large compared to the atomic scales); and finally we 
find the microscopic scale at MeV < A < tuq. 

On the other hand, assuming that the scale of i-colour is below the weak scale, bound 
states of the new sector are kinematically accessible to present and future collider experi- 
ments. However, since the SM particles are uncharged under i-colour, (quirk) loops would 
be required to couple both sectors leading to highly supressed production rates. Moreover, 
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from reference [2j , quirks are defined to be charged under (some) gauge groups of the SM 
(for example, charged under the electroweak interaction) quirks could be pair produced 
through electroweak processes. 

i — I 

Likewise, in Ref.[3] quirks are considered as vectors with respect to the electroweak 
gauge group without carrying QCD colour, but carrying i-colour charge. Therefore 
notice that there is no Yukawa coupling between the Higgs boson and quirks. Thus we 
discard the possible binding force which has been postulated for ultraheavy quarkonium 
taking over gluon-exchange. At this point, for the sake of clarity, it must be stressed that, 
taking simple assumptions, through this work we will focus on the case of uncoloured quirks. 

In this scenario, quirks can be copiously pair-produced at the LHC not through QCD 
interactions but via electromagnetic and electroweak interactions. As quirks would be 
long-lived particles as compared to the collider/detection time scale, different detection 
strategies can be undertaken according to the possible aforementioned micro, meso or 
macroscopic regimes. 

Finally, notice that possible quirky signals of folded supersymmetry in colliders have been 
studied in focusing on the scalar quirks (squirks). Contrary to usual supersymmetric 
partners of quarks, squirks (choosing a simple scenario) are expected to be uncolored, but 
instead charged under a new confining group, equivalent to i-QCD as introduced above. 
The study of spectroscopy performed in this paper is actually not sensitive to possible scalar 
nature of new fermions, and the main lines could be applicable to squirkonium as well. 

C. Exotic phenomenology 

According to [7], when the quirk pair is produced an excited bound state can be formed 
with invariant mass given approximately by the total center of momentum energy of the 
hard partonic scattering giving raise to the pair. In the microscopic regime this bound 
state would loose energy by emitting i-glueballs and bremsstrahlung towards low-lying 
states. Once they loose most of their kinetic energy these bound states (to be dubbed 
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Quirkonium) could decay via electroweak interaction. However, it is also possible other 
scenario in which the neutral and colourless quirk pair might have a prompt annihilation 
before they can loose energy enough to form the low lying quirkonium state. However, 
due to the non perturvative nature of the mechanisms it is difficult to estimate in advance 
the proportion of the events falling in each scenario, and therefore the possibility to 
detect these low lying states can not be discarded. Through this work, in particular we 
will focus on neutral bound states which can decay, e.g. to final-state dileptons, pro- 
viding a clean signature even admits a huge hadronic background as at the LHC experiments. 



New particles with a mass of up to several hundred GeV can be pair copiously produced 
at the LHC. One expects that quirks will be in general produced with kinetic energy quite 
larger than A. A significant fraction of this energy should be lost by emission of photons 
and i-glueballs prior to pair-annihilation. 

The two quirks will fly off back-to-back, developing a i-QCD string or flux tube. In usual 
QCD with light matter the string is broken up promptly by creating light quark- ant iquark 
pair; in i-QCD this mechanism is practically absent. The two heavy ends of the string would 
continue to move apart, eventually stopping once all the kinetic energy was stored in the 
string. The quirks would be then pulled together by the string beginning an oscillatory 
motion. 

Most examples of late-decaying particles that have been addressed in the literature yield 
missing energy, while quirks would annihilate into visible energy in most modes. Besides, 
as explained in 3 and A, only i-colour singlet states could be observed. 

As commented before, in a optimistic scenario, the excited bound state will emit 
i-glueballs and bremsstrahlung, towards low-lying states; then they would annihilate into a 
hard final state: di-lepton, di-jet, or di-photon. 



If we want to investigate whether or not it is possible to disentangle different state levels 
under the assumption of a given (large) quirk mass and a specific form of new i-colour 
interaction, then we could focus on a dimuon signal from the annihilation of a narrow 
resonance, since it is the most promising channel |7|; it then becomes crucial if the level 
spacing of different S'-wave states is enough to the foreseen mass resolution based on invariant 
mass reconstruction from a dimuon system. 
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According to the reference [7fl the detection of bound states at hadron colliders is reliable 
because the signal production is strong and peaked in invariant mass, and the dominant 
background is electroweak and diffuse. On the other hand, dimuon backgrounds from 
sources other than Drell-Yan can be suppressed by requiring no extra hard jets or missing 
energy. Besides, at the LHC, trigger and detector efficiencies are expected to be very high 
for high-mass dimuon events. 



Concerning other channels, it is expected that Quirkonium annihilation into a electron 
pair could be also a useful signal since the invariant mass peak is expected to yield a similar 
peak to dimuon but smaller and wider due to detector effects. 



On the other hand, in the squirkonium case, according to the reference jlOj the radiative 
decay (by soft radiation) from the highly excited states to the ground state can be ultimately 
detected by means of unclustered soft fotons in the uncolored case. Also, in ref. |9| it is 
pointed out the possibility to use the invariant mass peak of VT+photon, since this channel 
dominates the squirkonium annhilation at or near the ground state. If necessary, all these 
signals could aid to distinguish among different states. 



Focusing again on the dimuon signal, as it can be seen from ATLAS and CMS reports 
(llLQl one can consider a 2% accuracy for the transverse momentum of muons even at high 
momentum. Since the dimuon invariant mass should coincide to a good approximation, for 
small (pseudo)rapidities, with the transverse momentum error, 5(M W ) ~ 0.02 x Mqq 
(since p t ~ Mqq/2); letting Mqq vary along the range [100, 1000] GeV, the mass resolution 
should roughly take the values along the interval [2, 20] GeV. 

D. Model settings 

Hereafter we restrict our analysis to the range of A given by A << rriq but in the 
microscopic regime, namely 

few MeV < A < few GeV (1) 
As is well-known long ago, a non-relativistic treatment of the potential for conventional 
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heavy quarkonium has proved to be suitable on account of the asymptotic freedom of QCD. 
Moreover, one distinguishes between short- and long-distance dynamics of constituents in 
the bound state leading to an effective (static) potential of the type: 

V(r) = V s (r) + V L (r) (2) 

In particular we will write 

V(r)=-±+Br* (3) 

where the first term with fi = 1 would correspond to a Coulombic interaction, and the 
second one with v — 1 to a linear confining interaction. 

In this work we will consider firstly a Coulomb plus Linear potential (CpL) with [i = 1 
and v — 1; later we will use a more general Coulomb plus Power Law potential (CpP) 
with fi = 1 and v = 0.5, 1.5 as tentative possibilities. The motivation for the insertion of 
;hese power law potentials comes up from the clasical studies of Quarkonia (see Refs. |13|. 
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15].) in which are considered Coulomb like, linear and Cornell potentials, but as well 
the power law potentials are taken into account in order to cover possible deviations from 
them. In this way and focusing on the ease of Quirkonium, in reference fl a pure linear 
potential is taken into account, in reference [8] a Coulomb-like potential were considered. 
Therefore tracking the same philosophy than in the Quarkonia case possible deviations are 
also considered within this study. 

Moreover, the interaction accounting for the above static potential can be parameterized 
by the fermion (quirk) mass toq, where 100 GeV < uiq < TeV and an additional SU(Njc) 
gauge coupling (Njq stands for the i-colour number) which can be related with a i-colour 
scale A. 

In this work we will specify V(r) in Eq.Q as 

C a' 

V(r) = ar u - — (4) 

to be interpreted as CpL [y = 1) and CpP (v = 0.5, 1.5) potentials. Here a corresponds 
to the i-color string tension and we have introduced a i-color coupling a', alongside a 
group theory factor C, in close analogy to QCD potential models; hence, making a simple 
assumption, such group factor is taken as a mirror from QCD potentials (C = 4/3) and 
included into the infracolour coupling: i.e., in calculations we set oti co i our = Cot' . Of course, 
other numerical choices for C can be done, but a' depends on the A scale which is actually 



uncertain as we shall discuss later. 



On the other hand, in analogy to conventional QCD-inspired potential models, the i- 
colour string tension can be interpreted as a linear energy density (E/L), where E ~ A and 
L ~ A -1 . Hence the relation a ~ A 2 is expected to remain (approximately) valid, likewise 



the equivalent QCD expression a s ~ Aq CD (also derived from lattice calculations 16]), and 
finally 



a oc 



A ^ 



AqCD 



(5) 



In other words, a proportionality depending on the respective A and Aqcd between both 
string tensions could be expected from the above arguments. Basically, Eq.flH]) implies that a 
and A parameters are not independent of each other. For the sake of simplicity, the unknown 
proportionality factor will be set equal to unity, so that by fixing A one gets a (for given 
Aqcd and a s values). 

Focusing on the A scale, in this work in principle it corresponds to the microscopic 
A' scale depicted in Ref.jsj]. Numerically speaking, as it will be seen, the values of A 
were taken below and above of the QCD scale in a bandwidth; i.e. A = kAqcD with 
k = (0.1,0.4,1,10) to take into account the uncertainity about this quantity. In this 
way, the equation([5]) can be regarded as a comparison between the strength of the linear 
potential in both sectors SU(3)qcd and the new gauge group SU(3)i co i our , and it is 
intended to be an ansatz or a hint to determine numerically a proportionality between a s 
and a, in which subsequently the numerical uncertainity about the proportionality factor 
is diluted, taking into account the lack of knowledge about A. Besides, this comparison 
between diferent theories can be viewed to some extent in a similar way than in clas- 
sical physics, in which the strength of gravitational and the electrostatic forces are compared. 

Concerning the i-colour coupling constant, a', it would be related with A; as we are 
dealing with a non-Abelian i-colour binding force it implies a' is scale dependent. We will 
compute a' value at the running scale Q = 2mg according to [3|: 

/ / „x 12 7T 

a\Q) = (6) 



;il N IC - 2n Q ) In 



A- 
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where Njc is the i-colour number, and tiq the number of quirk generation at the running 
scale. From Eq.© and Eq.fjSJ), one can see that both parameters a and ol are depending 
on A, so that they are not independent quantities. Nonetheless, all those constraints have 
to be taken with a grain of salt and one should consider as well values deviating from those 
given in Eqs. as we will see later. 

In case of more quirk generations, additional active quirk should be taken into account at 
different energy scale thresholds. Nevertheless, for the sake of simplicity, and in view of still 
many unknowns in the different models, we will assume Nic = 3 and Uq — 1 throughout 
this work. 



As previously mentioned, we consider that the quirk mass lies in the range 100 GeV < 
rriQ < TeV. Therefore one can reasonably expect that the bound system indeed meets a 
truly non-relativistic regime, i.e. the relative quirk velocity v in the center of mass frame 
being substantially smaller than the value for bottomonium (v 2 ~ 10 _1 ). Focusing on 
quirkonium, a formal derivation of such non-relativistic limit from the relativistic degrees of 

fl n n 

freedom can be found in [8]; besides according to Reference [2] and 10], the bound state is 
formed in a highly excited state then it decays to the lower states loosing the main part of 
its kinetic energy. Therefore it is expected that in the lower levels near to the ground state 
this kinetic energy could be low enough to assume a non-relativistic aproximation. Also, as 
we shall see later, the numerical results obtained for the expected quirk velocity < v 2 > in 
the CoM frame justify this approximation. 



II. PROSPECTIVE SPECTROSCOPY OF EXOTIC STATES 

Since we are interested to perform spectroscopy for very heavy non-relativistic bound 
states, the Schrodinger radial equation must be solved: in a analytical way it could be done 
by means of a expansion of the quirkonium wave function in a complete basis; nevertheless 
here we will choose to solve it numerically, and therefore one should expect that the method 
to get the resulting mass spectroscopy followed in the QQ-onia package jl3| based on the 
resolution of the Schrodinger radial equation using the Numerov 0(h 6 ) technique, should 
work appropriately for our analysis of Quirkonium. However, we are confronted here to the 
lack of experimental data to set the ground state of quirkonium, in sharp contrast with, 
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e.g., the bottomonium or charmonium systems. Nevertheless, let us stress that in this work 
we are here mainly interested in estimating the mass spacing between different state levels 
rather than their absolute values. 

As commented in the Introduction, new interactions and particles can form very massive 
bound states. In this section we show the results for quirk (QQ) bound state system by 
sweeping through the scale range, few MeV < A < few GeV characterizing the i-colour 
force. 

First we will look at the results using a Coulomb plus Linear potential (CpL) (with 
fi = 1 and v = 1); later using a Coulomb plus Power Law potential (CpP) (with /i = 1 
and v = 0.5,1.5). Concerning the quirk mass, first we use vtlq = 100 GeV, and later 
rriQ = 500 GeV as representative values (nevertheless in some calculations we will take 
additional values). The energy levels, E n shown in Tables correspond to M n i = Iviiq + E n i 
where M n \ (or M) is the quirkonium mass level. In all cases, we set the ground state to 
be E\s = 0. Besides as relevant calculations we will display also the squared radial wave 
function at the origin (WFO) (or their derivatives), the size of each quirkonium level, and 
the mean value of the relative quirk velocity v in their center of mass frame, since it is 
used in some calculations besides, the obtained velocities will check the non-relativistic 
approximation. 
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A. CpL potential 



Let us start by considering the CpL potential. All Tables cited in this and successive 
sections can be found in Appendix. 

1. m Q = 100 GeV 



Results for the quirkonium spectrum with rriQ = 100 GeV and A = 0.1 Aqcd = 25 
MeV are shown in Table [H The corresponding parameters are a = 0.0018 GeV 2 and 
a'(Q = Zttiq) = 0.068. Concerning the mean radius, for the ground state we find a size 
similar to the Bohr radius ~ (m«') _1 , and increases for higher states as expected up 
to ~ ^fm; moreover, with these parameters we can found (85) states with sizes beyond 
2fm, which is in accordance with ~ A -1 . The quirk velocity in the CoM frame (v 2 ) ~ 10 -4 
slowly increasing with the n and I quantum numbers; these low v values plainly justify the 
non relativistic regime resulting from the QQ-onia package. 



We provide the squared WFO and derivatives divided by powers of the quirk mass ob- 
tained in our calculations, following the same behaviour with n and / as found in standard 
quarkonium (see for instance [17] and references therein). From the ground state WFO 



value we realize that IS state follows mainly a Coulombic behaviour |18|. This is not the 
case for higher resonances, for the P states case we find that a Coulombic (derivative) WFO 
underestimates the numerical value obtained from this potential. 

Let us stress that the energy level spacing (notably between S-wave resonances, of order 
of tens of MeV) would not permit the experimental discrimination by using the dimuon 
annihilation channel (and likely any other else). 

Table U shows the results for A = 0.4 • A QCD = 100 MeV (a = 0.029 GeV 2 and 
a'(Q = Ivtlq) = 0.08). The QQ level spacings turn out to be somewhat larger than in 
the previous case but still not enough to permit experimental discrimination. Something 
similar can be expected for A = Aqcd — 250 MeV as can be seen from Table IHIt with 
(a = 0.18 GeV 2 and a'(Q = Iraq) = 0.091). Here we find lower values for sizes of resonances 
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with respect to previous case (as expected since A increases). Besides, we observe a WFO 
value for the ground state somewhat greater than the expected for a Coulombic behaviour. 

The results shown in Table IIVI (appendix ) corresponds to the microscopic scale A = 
10 Aqcd- Here a'(Q = 2i71q) = 0.139, and a = 18 GeV 2 . The string tension turns to be 
much stronger than in the QCD case. The energy levels reach the GeV scale and the WFO 
grow to the ~ 10 3 GeV 3 values; according with previous trend the corresponding derivatives 
are growing also. The WFO value for the lowest state is ~ 3 times greater than the expected 
for a Coulombic behaviour. Concerning the level spacing this case is interesting since values 
among S'-wave states turns out to be of order of ~ 2 GeV, likely enough to be disentangled. 

2. m Q = 500 GeV 

We now set the quirk mass equal to 500 GeV, so quirkonium mass is of order of the 
TeV scale. In Table [V] the QQ spectrum is shown for A = Aqcd {p = 0.18 GeV 2 
and a'(Q = 2mg) = 0.073). In Table PVTl we show the results for A = 10 Aqcd- Here 
a'(Q = Ivtiq) = 0.102, and a = 18 GeV 2 . Again, as in the tuq = 100 GeV case at this 
scale, the level spacing among S'-wave states could be enough to distinguish experimentally 
these levels. 

B. CpP potential. 

Let us now give v values in the long-range term of Eq.fll]) different from unity. As in 
QCD, a larger (smaller) v leads to stronger (weaker) long-distance interaction. The general 
trends are similar to the CpL case seen in the previous section. 

1. ^ = 0.5 

Tables IVHI and IVIIII show the QQ spectrum for uiq = 100 GeV and tjiq = 500 GeV 
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respectively for A = A 



bottomonium case 



gcD 



As we can see the sizes of bound states are similar to the 



171 ]. 20]. Tables HXl and 1X1 (appendix) provides again the corresponding 
QQ spectrum for tuq = 100 GeV and mq = 500 GeV, but this time having set A = 
10 ■ Aqqd- Here we observe values of WFO for the ground state similar to the ones expected 
for a Coulombic behaviour; however higher resonances do not behave in this manner. 



2. u= 1.5 



To cover possible deviations from linear behaviour of the long distance part of the 
potential we analyze the CpP potential setting v = 1.5. Tables IXll and IXTTI (appendix) show 
the QQ spectrum for with A = Aqcd for rriQ = 100 GeV and vtlq = 500 GeV respectively. 
Tables IXIIII and IXIVI display the corresponding results for the QQ spectrum for itlq = 100 
GeV and m Q = 500 GeV with A = 10 Aqcd- 

In order to compare the effect of the above mentioned potentials, in Figure 1 we plot 
the nS level spacings A n s-is = M(nS) — M(1S) of quirkonium found with the CpL and 
CpP potentials (v = 1, 0.5, 1.5 respectively) for different itlq and A values. As far as we are 
interested in disentangling peaks of S'-wave resonances, it becomes apparent that this would 
be only possible in some cases (i.e. A = 10 Aqcd) where the level spacing is 0(1) GeV or 
larger. 

C. Other possible contributions from short distance potential 



Finally, to take into account other contributions which could be entangled in the short 
distance part of the potential, we consider higher (non perturbative) a' effective values. In 
order to consider this scenario, we do not use the Eq.([6]) for a' but we take it as a free 
parameter. On the other hand we keep the explicit dependence of A ( Eq.Q ) in a. In this 
case we also increase the A values from A = 10 • Aqcd up to A = 40 • Aqcd- 

By using QQ-onia code we find the results with itlq = 500 GeV for the mass level spacing 
A2S-1S = M(2S) — M(1S) which are shown in Table IXVI (Appendix). As we can see from 
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these situations, we find separation between levels tens of GeV, thus, in principle, we should 
be able to discriminate at least between these resonances. 



D. WFO vs. quirk masses. 

Next let us analyze the WFO dependence w.r.t. the quirk mass using the above explained 
CpP potentials. Here we focus on the IS ground level, by taking quirk mass values from 
rriQ = 100 GeV up to 500 GeV. Again we will take A = Aqcd, 10 Aqcd for each 
potential. For intermediate uiq values o does not change w.r.t. the mass values, However 
here a'(Q = 2mg) changes for each case according to Eq.Q. Figure 2 displays the obtained 
results. 

III. QUIRKONIUM DECAY 

Once computed the squared WFO, we can evaluate numerically the partial decay widths 
of neutral ( 3 Si) quirkonium ("0qq) to different final states. All of them are proportional to 
the ratio \Rs{0)\ 2 /M 2 , where M is the quirkonium mass. Subsequently we make estimates 
of the respective branching ratios (BR). 



We will follow a similar treatment as the authors of [3], [19] who considered the following 



where F( (Nic, Ri,eQ, SM) stands for functions containing the i-colour number, 
squared mass ratios Ri = Mf/M 2 (i = f,Z) , the quirk electric charge eg. The 
SM label means that those SM parameters involved in this calculation parameters are 
included. 




ipQQ decay modes: 




(7) 



• Decay to a 



pair 



W + W-) = F 2 (N IC ,Ri,SM) 



Rs(0)\ 2 



(8) 



M 2 
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Where F2(Nic,Ri,eQ,SM) stands for a function entangling the i-colour number, 
squared mass ratios Ri = Mf /M 2 (i = W, Z, tuq), and SM parameters. 



Decay to i-gluons ((/). Quirks couple to the i-gluon field of the SU(Njc) with coupling 
strength g' s = y/Arfa^, where a' s {Q 2 ) is given by Eq.([6]). 

n^ Q Q 9 99') = F 3 (N 2 C} c/ 3 )^^ (9) 
Here, F 3 (N 2 c ,a' 3 ) is a function of the i-colour number and the i-colour a' coupling. 

m>QQ -f 7 9 9) = F,(N IC ,ela' 2 ) l -^l (10) 



M 2 

\Rs(o)\ 



2 



T^ QQ — ► Z g'g') = F 5 (N IC ,ela' 2 ,SM)^j^ (11) 

To make the reading easy, the explicit form of F^i = 1, 5) coefficients can be found in 
references 

,0. 



A. Numerical results 



Once set the numerical values of parameters in the above expressions, the \R S (0)\ 2 values 
from Tables HI to IXIVI (appendix) allow one to compute the decay widths of ipQn(lS) to SM 
quarks (qq = uu, dd, ss, cc, bb, [ti if above the threshold]), leptons {11 = e ± , /i ± , r 1 * 1 ), 
and other boson decays {W^, 3g', ^y2g', Z2g'). The results are shown in Table IXVll for 
CpL and CpP (with v = 0.5, 1.5) potential using uiq = 100 GeV and itlq = 500 GeV at the 
above considered scales. 

In all cases the decay mode to SM quarks is the dominant channel. Decay to leptons 
shares roughly with a 33% for electron, muon and r pair respectively. As we can see, for 
A = Aqcd case if we take into account only the IS" decay, the total width is quite narrow 
~KeV, but we find similar values than in the heavy quarkonia case [21]. For A = IOAqcd 
case the total width increases roughly one order of magnitude. Nevertheless, if necessary, 



this analysis could be improved by adding upper nS levels contributions 22] : for instance if 
we compute the whole nS contribution using the m q = 500GeV^; A = IOAqcd case, with a 
CpL potential we find a total width pa 1.7 times the IS total width; using CpP v = 1.5 and 
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v = 0.5 potentials we find a factor ps 1.5 and ps 2. respectively. 



Concerning P-wave resonances (7^0), the corresponding widths satisfy 
so that, those contributions are suppressed with respect to the nS decays by a (D) factor 



D = 



1 l-R'j 



p 



2 



M 2 \R S (0)\ 2 

Taking n = 1 values in the CpL case from Tables , we find D ~ [10~ 5 , 10~ 4 ] in the 
toq = 100 GeV case and D ~ [10~ 6 , 10~ 5 ] for itlq = 500 GeV . Regarding the dependence of 
the BRs on the quirk mass (BR are independent of the ratio \R s (0)\ 2 /M 2 , but M enters also 
thorugh the functions Fj): in the range of interest 100 < tuq < 500GeV we find variations 
on the different BR less than a 1%. 

We can also check the BR variations with the i-colour number Njq'- by replacing for 
instance in the above expressions Njc = 6 -H- Njc = 3 -H- iV/c = 1, the .Bit! to bosons varies 
mainly 33% 17% 5% and the corresponding BR to SM quarks 55% <H> 68% ^ 78% 
(Pi? to leptons 12% <->■ 15% ^ 17% respectively). 



IV. CONCLUSIONS 

Spectroscopy of exotic states might play a fundamental role in the discovery strategy of 
new physics at the LHC and ILC In this paper we have focused on a simple extension of 
the SM, when a new SU (Nj) gauge group is added to the SM. The new interaction and new 
associated fermions have been dubbed i-color, quirks respectively. We assume that quirks 
are colorless, but otherwise carry SM quantum numbers, thereby coupling to gauge W ± , Z 
and 7 bosons. 

Quirks can bind forming very peculiar structures reminding. In this work we have focused 
on neutral QQ states called quirkonium, when the states are microscopic. We have performed 
a prospective study of quirkonium spectroscopy by employing a Coulomb plus Linear and 
Coulomb plus Power Law potentials as representative possibilities with parameters according 
to i-QCD requirements, as well as other effective contributions to analyze their impact. 
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Taking into account the wide range where the QQ bound state might be found, we 
have chosen the scale range, MeV < A < GeV with different i-colour A scales and quirk 
masses, finding sizes of several QQ n i resonances and their squared WFO values of states (or 
derivatives for / 7^ 0). We also extracted the level spacing among resonances using different 
A scenarios to determine whether or not it would be possible to discriminate different state 
levels. We also have computed total and partial decay widths. 
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Figure 1: Comparative plot of the nS energy levels w.r.t. the IS state, taken as the ground level, 
found with the CpL and CpP potentials [y = 1, 0.5, 1.5 respectively) for different run and A values: 
mlOO and m500 stands for the quirk mass value. L, 10L denotes A = (1, 10) Aqcd, respectively. 




Figure 2: Values of |i?is(0)| 2 /M 2 (in GeV) corresponding to the 15 level vs. M = 2toq (in GeV) 
using CpP potentials with v = 0.5; 1.5 (labeled as CPP05, CPP15 respectively). L, 10L denote 
A = (1, 10) Aqcd, respectively. The curve corresponding to v = 1 lies in between. 

V. APPENDIX: TABLES 



19 



Table I: Mass level spacings with respect to the ground state: A„;_i<j = M{nl) — M(1S) (MeV), 
using a Coulomb plus linear potential with rriQ = 100 GeV; A = 0.1 • Aqcd = 25 MeV; 
|i?2(0)| 2 /M( 2 + 2 ') (in GeV), and mean square radius (in fm). 



MeV. 



QQ LEVEL 


A„i_i5 (MeV) 

1 LI _L i_i \ / 


|p1?(o)I7m( 2 + 2 ') 


V(r 2 ) 

V \ ' 


(15) 





0.0092 


0.06 


(IP) 


153 


5.2 10~ 8 


0.23 


(25) 


154 


0.0012 


0.28 


(ID) 


183 


5.9 10~ 12 


0.45 


(2P) 


184 


2.1 10~ 8 


0.54 


(3S) 


185 


4.1 10~ 4 


0.58 


(AS) 


198 


2.2 10~ 4 


0.93 


(5S) 


206 


1.6 10" 4 


1.27 






« 10" 4 




II: The same 


as in Table fl] for a CpL potential with tuq = 100 GeV; 


A = 0.4 • Aqcd 


QQ LEVEL 


A nl _ ls (MeV) 


|/4(0)| 2 /M( 2 + 2 *) 


\/(r 2 ) 
v \ / 


(IS) 





0.016 


0.05 


(IP) 


229 


1.7 10" 7 


0.17 


(25) 


234 


0.0025 


0.21 


(ID) 


288 


8.3 lO -11 


0.28 


(2P) 


297 


1.1 10~ 7 


0.34 


(35) 


302 


0.0013 


0.37 


(45) 


344 


9.5 10" 4 


0.52 


(55) 


377 


7.9 10" 4 


0.65 



(v 2 ) « [10" 4 ,10- 3 ] 
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Table III: The same as in Table |U for a CpL potential with m Q = 100 GeV; A = A QCD = 250 MeV. 



QQ LEVEL 


A rd _ ls (MeV) 


\R l nl (0)\ 2 /M^ 2 + 21 ^ 




(15) 





0.025 


0.04 


(IP) 


348 


6.6 10~ 7 


0.12 


(25) 


375 


0.0059 


0.15 


(ID) 


473 


7.0 10- 10 


0.18 


(2P) 


508 


6.0 10- 7 


0.21 


(35) 


533 


0.0038 


0.24 


(45) 


649 


0.0031 


0.32 


(55) 


746 


0.0027 


0.38 



(v 2 ) w lO" 3 



Table IV: The same as Table H] for a CpL potential with m Q = 100 GeV; A = 10 • A QC d = 2.5 GeV. 



QQ LEVEL 


A nl _ ls (GeV) 


|P^(0)| 2 /M( 2 + 2 ') 




(15) 





0.218 


0.02 


(IP) 


2.51 


1.5 10" 4 


0.03 


(25) 


3.42 


0.129 


0.04 


(ID) 


4.15 


3.1 10~ 7 


0.045 


(2P) 


5.03 


2.1 10~ 4 


0.05 


(35) 


5.84 


0.109 


0.06 


(45) 


7.88 


0.099 


0.08 


(55) 


9.70 


0.092 


0.09 



(V 2 ) « lO" 2 
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Table V: The same as in Table H] for a CpL potential with m Q = 500 GeV; A = k QCD = 250 MeV. 



QQ LEVEL A ni _ 15 (MeV) |^(0)| 2 /M( 2 + 2 ^) 



(IS) 





0.058 


0.008 


(IP) 


831 


4.5 10~ 7 


0.04 


(25) 


850 


0.0063 


0.05 


(ID) 


1028 


2.0 lO" 11 


0.07 


(2P) 


1042 


2.3 10" 7 


0.09 


(35) 


1054 


0.0026 


0.1 


(45) 


1162 


0.0017 


0.15 


(55) 


1239 


0.0013 


0.19 



(v 2 ) « lO" 4 



Table VI: The same as in Table [J for a CpL potential with uiq = 500 GeV, A = 10 • Aqqd = 2.5 
GeV. 



QQ LEVEL 


A n/ _ 15 (GeV) 


|i4(0)| 2 /M( 2 + 2 ') 


7F) 


(IS) 





0.187 


0.007 


(IP) 


2.44 


1.0 10~ 5 


0.017 


(2S) 


2.98 


0.048 


0.022 


(ID) 


3.71 


2.0 10~ 9 


0.024 


(2P) 


4.18 


2.0 10~ 5 


0.029 


(35) 


4.64 


0.0350 


0.033 


(45) 


5.96 


0.030 


0.042 


(55) 


7.11 


0.027 


0.051 



(v 2 ) « 10- 3 
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Table VII: The same as in table U for Coulomb plus power law potential with v = 0.5, txiq = 100 
GeV; A = A QCD = 250 MeV. 



Drt T T?VTTT 
LJLJ Li-Cj V EjLi 


A-nl-lS (MeV J 


1 R z (Ci\\ 2 /A/f( 2 + 2 


1 l r 2\ 

v\ r / 


(15 J 


u 


n no7 


n n^ 
u.uo 


{if) 


/inn 

4UU 


y.o iu 


n 1 1 


[Zb) 


A 'XA 
4o4 


n nnf;s 


n i /i 

U. 14 




040 


r 7 1 n~ ii 


n 1 7 


(2P) 


579 


7.6 10~ 7 


0.21 


(35) 


604 


0.0039 


0.23 


(45) 


717 


0.0029 


0.32 


(55) 


803 


0.0023 


0.41 



(v 2 ) ps 10~ 3 



Table VIII: The same as in Tablefl]for a CpP potential with v = 0.5, itiq = 500 GeV; A = Aqcd = 
250 MeV. 



QQ LEVEL 


A ni „ 15 (MeV) 


\R l nl (0)\ 2 /M( 2 + 21 ) 




(15) 





0.060 


0.01 


(IP) 


880 


5.6 10~ 7 


0.04 


(25) 


911 


0.0072 


0.05 


(ID) 


1112 


6.0 10~ 12 


0.07 


(2P) 


1132 


3.0 10~ 7 


0.08 


(35) 


1150 


0.0031 


0.09 


(45) 


1276 


0.0020 


0.14 


(55) 


1363 


0.0015 


0.18 



(v 2 ) « 10- 4 
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Table IX: The same as in Tablefl]for a CpP potential with v = 0.5, uiq = 100 GeV; A = 10-Aqc^ = 
2.5 GeV. 



QQ LEVEL 


A-, i9 (GeV) 

III — ±lJ \ / 


|^,(0)| 2 /M( 2 + 2 ') 


J(r 2 ) 


(15) 





0.539 


0.014 


(IP) 


5.09 


6.3 10" 4 


0.024 


(25) 


6.94 


0.298 


0.031 


(ID) 


8.36 


3.1 10" 6 


0.032 


(2P) 


9.84 


0.0013 


0.040 


(35) 


11.28 


0.227 


0.046 


(45) 


14.63 


0.190 


0.060 


(55) 


17.42 


0.167 


0.073 



(V 2 ) « 10- 2 



Table X: The same as in Table[I]for a CpP potential with v = 0.5, rriQ = 500 GeV; A = 10- Aqcd = 
2.5 GeV. 



QQ LEVEL 


A n/ _ 15 (GeV) 


|P^(0)| 2 /M( 2 + 2 ') 


7F> 


(15) 





0.298 


0.007 


(IP) 


4.60 


9.0 10~ 5 


0.012 


(25) 


6.14 


0.107 


0.015 


(ID) 


7.39 


5.0 10" 8 


0.016 


(2P) 


8.44 


1.0 10" 4 


0.020 


(35) 


9.61 


0.075 


0.023 


(45) 


12.21 


0.061 


0.031 


(55) 


14.33 


0.052 


0.038 



{v 2 ) « 10- 3 
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Table XI: The same as in Table U for a CpP potential with v = 1.5, mq = 100 GeV; A = Aqcd = 
250 MeV. 



QQ LEVEL 


A n , io (MeV) 

fit — IO \ * / 


\R l ,(0)| 2 /M( 2+2 


\Z(r 2 ) 


(15) 





0.022 


0.05 


(IP) 


313 


4.6 10~ 7 


0.13 


(25) 


329 


0.0049 


0.16 


(ID) 


416 


1.5 10~ n 


0.20 


(2P) 


445 


4.4 10~ 7 


0.23 


(35) 


464 


0.0034 


0.25 


(45) 


568 


0.0029 


0.33 


(55) 


662 


0.0027 


0.39 



< V 2 > « IO" 3 



Table XII: The same as in Table [I] for a CpP potential v = 1.5, mg = 500 GeV; A = Aqcd = 250 
MeV. 



QQ LEVEL 


A nZ _ 15 (MeV) 


|P^(0)| 2 /M( 2+2 




(15) 





0.058 


0.01 


(IP) 


813 


3.7 10~ 7 


0.04 


(25) 


827 


0.0058 


0.05 


(ID) 


991 


2.0 10~ 12 


0.08 


(2P) 


998 


1.8 10~ 7 


0.10 


(35) 


1005 


0.0022 


0.11 


(45) 


1094 


0.0015 


0.16 


(55) 


1198 


0.0011 


0.20 



{v 2 ) « 10- 4 
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Table XIII: The same as in Table fl] for a CpP potential with v = 1.5, tuq = 100 GeV; A = 
10 • Aqcd = 2.5 GeV. 



(cj/l,/ LiHj V CjLi 


A-nl-lS l^ eV J 


1 R z (TlM 2 /A/f( 2 + 2 


//„2\ 

v\ r / 


{Lb) 


u 


nils 


n D97 




119 
1.1Z 


/ill n~ 5 

4.1 1U 


n n/i i 


f9 


1 88 


U.UD f 


U.UO / 




9 "il 


o. i iu 


n nfin 

U.UDU 


(2P) 


2.83 


5.8 10~ 5 


0.070 


(35) 


3.27 


0.059 


0.079 


(45) 


4.52 


0.056 


0.097 


(55) 


5.69 


0.054 


0.112 



(v 2 ) « 10~ 2 



Table XIV: The same as in Table fl] for a CpP potential with v = 1.5, rriQ = 500 GeV; A = 
10 • Aqcd = 2.5 GeV. 



QQ LEVEL 


A ni _ 15 (GeV) 


|p4/(o)I7^ (2+2 ° 




(15) 





0.160 


0.008 


(IP) 


1.73 


5.0 10~ 6 


0.023 


(25) 


1.91 


0.026 


0.029 


(ID) 


2.04 


8.0 10~ 10 


0.031 


(2P) 


2.62 


5.0 10~ 6 


0.040 


(35) 


2.80 


0.019 


0.044 


(45) 


3.51 


0.016 


0.057 


(55) 


4.15 


0.015 


0.067 



{v 2 } ps 10" 3 
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Table XV: Mass level spacing A 2S ~is = M(2S) - M(IS) (with m Q = 500 GeV) using A = 
(10,20,40) • A.qcd and different a' values. 





A 


a' 


A 25 _ 15 (GeV) 


10 


' Aqcd 


0.3 


« 10 


20 


• Aqcd 


0.3 


» 12 


40 




0.6 


» 30 



Table XVI: Branching Ratios of iPqq{1S) (M = 2m Q = 200 GeF) and ^qq{IS) (M = 1000 GeF) 
to SM quarks (gg),leptons (ZZ), and other boson decays. r(CpL,CpP15,CpP05) stands for the total 
decay width (in KeV) using each potential. 



i> QQ (lS) (M = 200 GeV) 


A(MeV) 


r(CpL,CpP15,C P P05)(KeT/) BR(ll)(%) BR(qq){%) BR(Boson)(%) 


25 


(31, , ) 17 76 7 


100 


(54, , ) 16 76 8 


250 


(86, 78, 94) 16 75 9 


2500 


(833,451,2060) 15 68 17 


4>qq(1S) (M = 1000 GeV) 


A(MeV) 


r(CpL,CpP15,C P P05)(KeF) BR(ll)(%) BR(qq){%) BR(Boson){%) 


250 


(165,164,169) 16 80 4 


2500 


(555,474,882) 16 76 8 
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